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Abstract. Fluid reductions of the Vlasov-Ampere equations that preserve the 
Hamiltonian structure of the parent kinetic model are investigated. Hamiltonian 
closures using the first four moments of the Vlasov distribution are obtained, and 
all closures provided by a dimensional analysis procedure for satisfying the Jacobi 
identity are identified. Two Hamiltonian models emerge, for which the explicit closures 
are given, along with their Poisson brackets and Casimir invariants. 
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1. Introduction 


The Vlasov-Ampere set of equations is a suitable framework for describing the dynamics 
of systems interacting through electrostatic forces. In this work, we focus on the study 
of electrostatic plasmas even though the results may be applied to more general systems 
described in part by the Vlasov equation. We consider a one-dimensional plasma made 
of electrons of unit mass and negative unit electric charge, evolving in a neutralizing 
background of static ions. The evolution of the distribution function of the electrons 
/, defined on phase space with coordinates (x,v), and electric field E is given by the 
Vlasov-Ampere equations, 


( 1 ) 

( 2 ) 


d t f = -vd x f + Ed v f, 
d t E = —47rj, 


where E and j are the fluctuating parts of the electric field E and the current density 


j = — J vf dv respectively. We assume vanishing boundary conditions at infinity in the 


velocity v so that integrals such as the charge and current densities are well-defined. In 
this work, we limit ourselves to the study of systems of unit length in the spatial domain 
x with periodic boundary conditions. The fluctuating part of the electric field is defined 


by E = E — fy E dx. The system is fully nonlinear, but has a form that builds in the 


preservation of the spatial average of E and maintains momentum conservation. 

The use of fluid reductions to describe the dynamics of a plasma is ubiquitous 
in plasma physics. Indeed, this usually allows one to decrease the complexity of the 
problem at hand and to gain physical insight into the phenomenon under investigation 
since the dimension of phase space is reduced. Fluid reductions of the Vlasov-Ampere 
equations are done by introducing fluid quantities such as the fluid moments 



( 3 ) 


The associated dynamical equations are then obtained by multiplying Eq. (P by v n and 
integrating with respect to the velocity. This leads to 


( 4 ) 

( 5 ) 


d t P n = -d x P n+ 1 - nP n _iE, 
d t E =4vr P\, 


't r n ~ 


for all 7i G N. In order for this system to be reduced, one has to truncate the 
infinite sequence of Eq. (14|) . Truncating this system at order N , that is considering 
{Po, Pi, -, Pn, E) as dynamical field variables, one can see from Eq. ([4]) that the time 
evolution of Pn depends on Pn+i- As a consequence, it is necessary to express Pn+i in 
terms of (P 0 , Pi,..., Pjv, E) in order to close Eqs. (J4J and (JSJ) and thereby obtain a fluid 
reduction. 

Many models have been proposed based on as many closures with various 
requirements (see, e.g., Refs. HI 121 Eli). A usual procedure consists in assuming a 
particular form for the distribution function / (e.g., Dirac, Maxwellian,...) depending 
on a finite number of parameters, and expressing the closure with respect to these 
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parameters [5]. Alternatively, closures have been constructed in order to recover certain 
kinetic effects 0 0 0 ej noi m Eg. In any event, a reduction by closure should 
be such that, if the parent model possesses a Hamiltonian structure PH EH H51 HE], 
then the resulting fluid model should also have one, after discarding all the terms 
that are supposed to provide dissipation. A closure procedure ignoring this aspect 
could potentially lead to the introduction of some nonphysical dissipation nzi m- 
Consequently, here we use a procedure that preserves the Hamiltonian structure of the 
parent kinetic (Vlasov-Ampere) system, which is one of its most important structural 
features. Specifically, in this work we present a model for the first four fluid moments of 
the distribution function, namely the density p, the fluid velocity u, the pressure P and 
the heat flux q. This allows us to account for the time evolution of the heat flux, which 
is of great importance for the study of transport phenomena inside the plasma. For such 
a model with four moments, one has to find a closure for the fifth order moment of the 
distribution function, namely P 4 . Here, we determine all the closures, obtained from 
a procedure based on dimensional analysis, that preserve the Hamiltonian structure of 
the parent model pun given by Eqs. ([Tj) and ([2]). We show that there are only two 
such Hamiltonian closures. The equations of motion of one of these two models are 
identical to the ones obtained with a bi-delta reduction [20j [21, 22], i.e., assuming that 
the Vlasov distribution has the form 

f(x,v,t) = uji5(v - nf) + uj 2 S(v - p 2 ), 

where a+ 2 and p ij2 depend on space and time. It should be noted here that we obtain 
these equations without any assumption on the special form of the distribution function. 
We provide the explicit expressions of the Hamiltonian and the Poisson bracket for the 
two Hamiltonian models. In addition, we derive the global Casimir invariants, which are 
specific invariants resulting from the knowledge of the Poisson bracket. These conserved 
quantities can be used, e.g., to ensure the validity of a numerical simulation of the 
equations of motion. 

The paper is organized as follows. In Sec. [2] we describe the methodology used for 
the derivation of the two Hamiltonian reduced models. We start from the definitions 
of the appropriate variables, namely, the reduced fluid moments. Subsequently, we 
introduce our method, based on dimensional analysis, which leads to models that obey 
the Jacobi identity. We show that there are only two such models. In Sec. [3l we 
analyze the two resulting Hamiltonian closures, providing explicit expressions for their 
Hamiltonians, Poisson brackets, and Casimir invariants. 

2. Method 

2.1. Reduced moments 

Our purpose is to build a Hamiltonian fluid model for the first four moments of the 
distribution function, namely the density p, the fluid velocity u, the pressure P, the heat 
flux q and the electric held E. These models will be referred to as 4+1 held models, 
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where the 4 refers to the four first moments of the Vlasov distribution (or equivalently 
to p, u, P and q ) and the 1 refers to the electric held E. We begin by considering the 
Poisson structure of the parent model with (/, E) as dynamical held variables. It was 
shown in Ref. [23] that the system of Eqs. o-o possesses a Hamiltonian structure 
with Poisson bracket 


{F,G} 


f 


d x F f d v Gf - d x G f d v F f + 4t r(F E d v G f 


G E d v F f ) 


dxdu, 


( 6 ) 


where Ff (resp. F E ) denotes the functional derivative of F with respect to / (resp. E). 
In addition, Bracket (|6]) is bilinear and satishes the Leibniz rule and the Jacobi identity. 
The Hamiltonian of the system is given by 
f v 2 f E 2 

H = j f — dxdv + J — dx, (7) 

where the hrst term accounts for the kinetic energy of the particles and the second 
one corresponds to the energy of the electric held. Together with Bracket (]6|) , this 
Hamiltonian leads to Eqs. ([Q) and (J2]) by using d t f = {/, H} and d t E = {E.H,}. 
We recall that such a bracket has Casimir invariants, i.e., functionals C that Poisson- 
commute with any other functionals of the Poisson algebra, {C, F} = 0 for all F. 
Bracket (J6]) has the following global (i.e., independent of the coordinates x and v ) 
Casimir invariants 

Ci = J <p(f) dxdv , 

C’2 = [ E dx, 


for any scalar function <p, and a local Casimir invariant 


C L = 


d x E + 47t 



which is equivalent to Gauss’s law. 

The change from the kinetic to the fluid description is done by performing the 
change of variables defined by Eq. (J3]) in Bracket ([ 6 ]) and Hamiltonian ([TJ) . The latter 
becomes 



Making use of the chain rule to transform the functional derivatives, Bracket ([ 6 ]) becomes 

mmm 


{F,G} 



Pi+j-i(Gjd x Fi 


Fjd x Gi) + 4tt Pj-^GjFE - F 3 G E ) 


dx, 


( 8 ) 


where F n denotes the functional derivative of F with respect to P n , and summation is 
implicit over the repeated indices i and j. Because we want to construct a Hamiltonian 
model for the first four moments of the distribution function, we consider functionals of 
the kind F[P 0 , P 1 , P 2 , P 3 , E]. However, the Poisson bracket (jBJ) of two functionals of this 
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kind depends explicitly on two additional moments, namely P 4 and P 5 . In order to close 
the system, these two additional moments need to be expressed in terms of P n < 3 and 
E. As a result, the Jacobi identity is no longer satisfied in general, and the resulting 
truncated and closed bracket is not of Poisson type. Consequently, the resulting system 
is not Hamiltonian, or in other terms, the reduction procedure potentially includes 
dissipation. We notice that the closure has to be performed on two moments, P 4 and 
P 5 , which slightly differs from what has been stated in the introduction, concerning 
the closure performed on the equations of motion directly, where only one additional 
moment, P 4 , needs to be closed. However we shall see in Sec. 12.21 that the expression of 
P 5 is entirely determined by P 4 . 

We introduce the reduced fluid moments, which we find to be more suitable variables 
for our purpose, 

P = f f u = ~ p / v f S n = J {y- u) n f du, (9) 

for all n > 2. The first and second ones correspond respectively to the usual density and 
fluid velocity. The higher-order moments are the central fluid moments with a specific 
scaling with respect to the density. The change from the usual fluid moments P n to the 
reduced fluid moments ( p,u,S n ), used hereafter, is invertible so that the results, even 
though they are expressed in a different set of coordinates, are equivalent. This change 
is given by 



for all n > 2. The inverse of this transformation is given by 


P 0 = p, Pi = pu , P n = P 


u 


+ £ 

m =2 


n 


m 


p m u n ~ m S rr 


Explicitly for the first four moments of the distribution function, this change of variables 
is given by 

Pi 


P = Ph 


u = 


Pn’ 


S3 = T (A - 3A? + 2 A ], 


p§ 


pi 


Pn 


' p2 

M) 


with the inverse 


P 0 = p, Pi = pu, P 2 = p ( u 2 + p 2 S 2 ) , 


P3 = P ( u 3 + 3p 2 uS 2 + p 3 S 3 ) . 


In terms of the moments, Hamiltonian ([71) is 

n = \f ( 'p u2 + p 3s2 + dx - ( 10 ) 

The first part of Hamiltonian (TTOl) accounts for the kinetic energy of the system while 
its second part corresponds to the internal energy. The last term, which accounts for 
the electric energy, remains unchanged compared to Eq. ([71) . By considering functionals 
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of the kind F\p, u, S 2 , S 3 , E} and using the chain rule for the functional derivatives (see 
Appendix C for more details), Bracket 


{F, G} = 


G u d x F p 


takes the form 
F u d x G p + 4n(G u F E -F u G E ) 
FiG 


- (G u Fi — F u Gi)d x Si + a.ij -- + d x f— ^ /3 *j —- da:, (11) 

P P P \P J P J 

where F) denotes the functional derivative of F with respect to S). From now on and 

unless otherwise stated, summation from 2 to 3 over repeated indices is implicit. The 

matrices a and ft have indices ranging from 2 to 3 such that 


a = d x 


2 S 3 2 S 4 - 3S% \ 

3S 4 - 6S% 3S 5 - 12S' 2 S , 3 ) ’ 


P 


AS 3 5S 4 — 9Sf ) 

5S 4 - 9 Si 6S 5 - 24 S 2 S 3 J J 


We notice that d x /3 = cc + cd, a property that ensures that Bracket (HI]) is antisymmetric. 
From Definitions (fT21) we see that the closure requires reexpression of S 4 and S 5 , i.e., 
one has to express these two reduced moments with respect to the dynamical variables 
(p,u,S 2 ,S 3 ,E) such that Bracket (HID satisfies the Jacobi identity. 

We remark that Bracket (fill) has several subalgebras. Trivial ones include F[p] (i.e., 
the algebra of functionals of the type F[p\), F[u], F[E ], F[p, E], and non-trivial ones 
include F\p,u], F[p, S 2 , S 3 ], F[u,E], F[p,u, S 2 , S 3 ], F[p,u,E} and F\p, S 2 , S 3 , E\. The 
most interesting one is the subalgebra of functionals F\p,S 2 ,S 3 ] for which p becomes 
a Casimir invariant. The existence of this subalgebra is the reason for considering the 
reduced fluid moments S n . 


2.2. The Hamiltonian constraints 


In order to be a Poisson bracket, Bracket (ITTT) must satisfy the Jacobi identity, 

{F, {G, H}} + {H, {F, G}} + {G, {H, F}} = 0. 

Here we determine the conditions on S 4 and S 5 resulting from requiring the Jacobi 
identity. We begin by assuming that S 4 and S$ depend on p, u, S 2 , S 3 , E and their 
derivatives dfp, d x u, d x S 2 , dfS 3 , d x E for n lower than some order v. Using the result 


obtained in Appendix A[ we conclude that S 4 and S 5 do not depend on p, u, E and 
their derivatives d x p, d%u, dfE. In addition, we show in Appendix B that in order for 
the Jacobi identity to be satisfied, we need to impose 

r iljmTfkij ^ikjnrdlij t 

for all z, fc, / and m ranging from 2 to 3, where the summation is implicit on j, and 

For instance, forZ = 2,m = 3,i = 2 and k — 3, we end up with 72337323 = 0 since 
7222 = 0 and 7223 = 0 for u > 2. From Eq. 0 . we have 37233 = 27323 , therefore 
7233 = 0, or equivalently 
dS 4 


dd v - l S, 


= 0 . 
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Using Eq. (IB. 51) leads to 
da 23 


S 3 


= 0. 


ddfS 2 

Since this has to be true for any value of S 3 , we thus conclude that 7232 = 0, i.e., 

dS 4 


dd"- ] S 2 


= 0. 


Concerning S 5 , Eq. (IB.5|) for l — i — 3 leads to 

dS 5 


fikj 




= 0. 


(13) 


There are two solutions to Eq. (fT3l) . The first solution is given by 
dS 5 dS 5 


dd^S, 


= 0, 


>2 dd^S 3 

The second solution requires det (3 = 0, which, using Eq. (1I2)) . can be written as 

(5S 4 - 9S 2 ) 2 


= 0. 


S 5 = 4S 2 S 3 + 


4S 3 


Since S 4 does not depend on df. 1 S 2 and df. l S 3 , we again have 
dS§ dS§ 


= 0, 


= 0. 


ddf-hS-i ddf- l S :i 
In what follows we will see that the second solution does not lead to a Hamiltonian 
closure. By induction on v down to v = 2 we show that S 4 and S 5 have to be functions 
of S 2 and S 3 only. These conditions are necessary but not sufficient, i.e., for any functions 
S 4 and S 5 of S 2 and S 3 , Bracket (fill) does not satisfy the Jacobi identity in general. 

We compute in Appendix C the necessary and sufficient conditions on the closures 
for a fluid bracket of the type (HIT) to satisfy the Jacobi identity. For four fluid moments, 
Eqs. (1C.91) and (1C. 101) are 

'<9S 4 


dS 5 _ dSi 

— + 


dS 2 
dS 5 
dS 3 


dS 4 

dS, 


dS 3 

2 

+ 


6 S 5 = 4S 3 3S 2 + 


dS 2 

dS 4 

dSf 

dS 4 

dS~ 2 


3 Sc 


dS 4 


(9S 2 2 - 5S 4 ) 


(14) 

(15) 

(16) 


We see from Eq. (UB1) that the expression of S 5 is fully determined by S 4 . By introducing 
the expression for S 5 given by Eq. (fl 6 |) into Eqs. (1141) and (fT5l) . we end up with the 
following two nonlinear second order partial differential equations: 

<9 2 S 4 


4S; 


4S a 


d 2 S A 

~dSf 

d 2 S 4 

dS 3 dSo 


dS 2 dS 3 

d 2 S 4 


(9S 2 2 - 5S 4 ) - d f n 4 d f n 4 = 12 S 3 , 


dS 2 dS 3 


(9Sf - 5S 4 ) + 12S 2 = ( 


dS 4 

dS 3 


+ 2 


ds 4 

dS 2 


(17) 


(18) 


Provided that these two equations are satisfied, Bracket (HI]) is a Poisson bracket and 
the resulting system is Hamiltonian. Solving these equations in general is challenging; 
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consequently, in what follows we restrict ourselves to the set of solutions provided by 
dimensional analysis |27j . 

2.3. Closures based on dimensional analysis 

We consider all the closures for the fifth-order moment S 4 = g(S 2 , S 3 ) that satisfy the 
constraints given by Eqs. (1171) and (lT 8 l) based on a dimensional analysis argument. In 
order to proceed, we assume that the closure S 4 = g(S 2 ,Sz ) does not depend on any 
further dimensional parameters. This would not be the case for, e.g., diffusion-like 
closures (Fourier’s law, Fick’s law, etc...) that introduce phenomenological parameters 
resulting from various hypotheses based on characteristic scales of the dynamics of the 
system. Indeed, in diffusion processes, diffusion coefficients replace information on the 
particle interactions, thus removing small scale dynamics. Instead, we would like our 
reduction procedure to be very general and not to depend on the geometry of the system. 
Consequently, we seek Hamiltonian closures where S± = g(S 2 , S 3 ) do not depend on any 
further dimensional parameters. 

It can be shown from Eq. ([9]) that the dimensions of the Sfl s, denoted [A n ], are not 

independent. Indeed, for all n > 2 we have [S n ] = A n , where A = L 2 T ^ 1 with L and 

T denoting the units of length and time, respectively. As a consequence, the closure 
S4 = g(S2, S3) involves three quantities and a unique physical dimension A. Making use 
of the Buckingham 7 r theorem [27], there exists two dimensionless quantities, denoted 
( and £, such that S 4 = g(S 2 ,S 3 ) reduces to ( = _R(£). Therefore, this procedure 
eliminates one of the variables in the closure. Defining ( = S4/S 2 and £ = S3/S ^' 2 and 
inserting these expressions into Eqs. flT7|) - (fT8]) . we get the following two constraints: 

3£A"(6£ 2 + 9 - 5R) + R'(3£R' - 18£ 2 + R ~ 9) + 16£A = 24£, (19) 

R'\ 6£ 2 + 9-5 R) + R\R' - 5£) +4 R = 12. (20) 

To solve Eqs. (fT9l) - (120l) . we compute their values for £ = 0. Defining Ro = R( 0), 
R ' 0 = R'( 0), and Rq = R"(0), Eqs. (1T9|) and (120]) become 

R' 0 (R 0 - 9) = 0, (21) 

R"(9 - 5R 0 ) + R ' 0 2 + 4R 0 = 12. (22) 

Equation (1211) has two solutions: R ' 0 = 0 and Ro = 9. Equation (1221) then reads 

Rq = 4(3 — Rq)/( 9 — 5R 0 ) for R ' 0 = 0 and R.f = 12(3i?g — 2) for R 0 = 9. We now 

differentiate Eqs. (fT9l) and (j20l) with respect to £ and evaluate them at £ = 0. This gives 
us 

R"( 9 - 7Rq) + 2Rq + 8f? 0 = 12, (23) 

R'o (9 - 5R 0 ) - R' 0 (3Rq + 1) = 0. (24) 

Using R 0 = 9, Eq. (1231) together with Eq. (1221) leads to Rq = —2/3 and R'f = —48. 

As a consequence, this solution is not real and of no interest for our purpose. The 
other solutions satisfy (R 0 , R' 0 , Rq, Rq) = ( 0 , 0 ,4/3, 0 ) and (R 0 , R' 0 , Rq, Rq) = ( 1 , 0 , 2 , 0 ), 
where Rq = R'"( 0). Since the solution is unique for a given set of initial conditions, 
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there exist only two solutions to Eqs. (1191) and (19UD . Moreover, one can see that R(—£) 
is also a solution of these equations. Consequently, the two solutions R are even with 
respect to £. These two solutions are described in the next section. 


3. Hamiltonian fluid models with 4+1 fields 


3.1. Model without normal variables 

We consider the first solution to Eqs. (fT9j) - (l20j) . corresponding to (Rq, R! 0 , R'f Rff) = 
(0,0,4/3,0). As mentioned in Sec. 12.31 the solution R is even. Thus we introduce 
R(£) = R(r /), where 7] = £ 2 . Then, Eqs. (1T9]) - (120]) become 

f [3r]R"(6ri + 9-5 R) + R'(9 + 3 rjR' - 7R) + 4R - 6] = 0, (25) 

2riR"(6ri + 9-5 R) + R \9 + 2r)R’ - 5R + rj) + 2R - 6 = 0. (26) 

By linearly combining these equations to eliminate terms in R", and introducing the 

new variable fi — —(R — 3?/ — 9)/5, we end up with an Abel equation of the second kind 

(see, for instance, Ref. [28]): 


fi/a — /a = 


67] + 24 
25 ’ 

which has the parametric solution 

. (2 — 5r) 2 


t](t) = K 


(3 - 5r) 3 


4, ]i{r) = Kt 


(2 — 5 r) 2 


(27) 


(3 — 5r) 3 ’ 

where K is some constant to be determined. Inserting Eq. (]27l) into Eqs. fl25 ]) -([26 p 
implies these equations are satisfied if and only if K = 27. This leads to an explicit 
expression for the closure ( = S±/S% = R(£) given by 

n [c| O I 4 + 4f (0 2 - - 8) - sap + - m? - 2) - 4g] ^ (28) 


where 


m = 


[i-2^+3 m+m 2 f 

\/e^+e) 3 -2-ioe+e 


1/3 


Furthermore, by using Eq. (TT6T) . S 5 is given by S 5 = S 2 S 3 T(f) with 

- 2 flK) _3e-9 ■ (29) 

In summary, the Hamiltonian and the Poisson bracket resulting from this closure are 
given respectively by Eq. (11(71) and Eq. dill) with a and f given by Eq. (1191) with 

S 3 


s 4 = S z 2 R 


sy 
s . 


(30) 


s: 


3/2 


S 5 = s 2 s 3 t 


(31) 
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Figure 1. Color map of S 4 (in A 4 units as defined in Sec. 12.311 given by Eq. (l30l) as a 
function of S 2 (in A 2 units) and S 3 (in A 3 units). 


S 3 


4 

2 

0 

-2 

-4 



xlCT 2 

8 

4 

0 

-4 

-8 


1 


2 3 4 5 

S 2 


Figure 2. Color map of S 3 (in A 5 units as defined in Sec. 12.31) given by Eq. m as a 
function of S 2 (in A 2 units) and S 3 (in A 3 units). 


where R and T are given by Eqs. (|28h - f[29T) . The dependence of the functions R and 
T in their arguments is not trivial. In order to help the reader visualize the closure 
relations corresponding to Eqs. (l30]l - (l3Tf . we provide, in Fig. |T] (resp. Fig. [ 2 ]), color 
maps showing the dependence of 64 (resp. S 5 ) on S 2 and S 3 . As a side note, we 
remark that as S 3 goes toward 0, 64 also goes to 0, as shown in Fig. Q] Thus, with this 
closure, non-asymmetric distribution functions cannot exist. Consequently, the physical 
relevance of this solution is questionable. 

In order to further characterize this Poisson bracket, we investigate its Casimir 
invariants. These are functionals C[p,u, S 2 , S 3 , E] that commute with any other 
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functionals F, i.e., {F,C} = 0 for all F. In particular, C commutes with the field 
variables. As a consequence, we must impose 

{p, C} = -d x C u = 0, 

which leads to 

C[p, u, S 2 , S :i , E] = Ki J u dx + D[p, S 2 , S 3 , E ], 
where K\ is constant. Imposing that C commutes with u leads to 

{u, C} = —d x D p - 47 tD e + -pDid x Si = 0, 
whose solution is given by 

D\p, S 2 , S 3 , E\ = K 2 J p<f(S 2 , S 3 ) dx + K 3 J E dx, 

where K 2 and K 3 are constant. Imposing that C commutes with S t leads to 

{Si, C} = ~—d x Si + K^cfj - — ifijcfj) = 0, (32) 

p p p 

where <pi = d(f/dSi. We then solve the associated homogeneous equation (A\ = 0). 
Again making use of the Buckingham 7r theorem, we assume that there exist a real 
number a and a function if such that (f = S%if(S 3 /S 2 ^ 2 ). The resulting equations are 

8(a - l)a^if + if' (3 - 18a + 30£ 2 - 24a£ 2 -3 R+ 10 aR - 9 fR') 

+ 3£ (9 + 6£ 2 - 5 R) if" = 0, 

2 aif + if'(-9£ + 4a£ + 3 R') + (-9 - 6£ 2 + 5 R)if" = 0, 

4a^[3 - 9a + (-1 + 5 a)R - 3fR'] = 3f [^'(1 - 4a + (-17 + 20a)A 
- 8(a - 1 )T - 6 fR' + 6 £T') - 3£(7 + 5 R- 4 T)if"}, 

AafR 1 + if'[3 — 18a + 5(—3 + 2 a)R + 6 T — 6 fR' + 6 fT'] 

= 3£(7 + 5 R- 4T)if". 

Combining the first two equations leads to 

2a(4a - l)£if + [3(1 - 6a + £ 2 - 4 a£ 2 ) + (-3 + 10 a)R]if = 0, 
whose solution is given by 

. , . , ( f 2ay(4a — 1) \ 

m = V-oexp I j 3(6o _ 1} + mia _ 1) + R{y){3 _ 10a ) ^ I ’ 

where ifo is a constant. Inserting this expression into the previous equations provides the 
necessary constraints a = 0 and if(f) — ifo- As a consequence, this model does not have 
Casimir invariants of the entropy-type [29J, i.e., of the form J pf>(S 2 , S 3 ) dx. Moreover, 
we can show in a similar way that Eq. (l32l) has no solution for the nonhomogeneous 
case (Ii 1 0), i.e., K\ = 0 is required. The Poisson bracket resulting from this closure 
has only two global Casimir invariants, given by 

C\ = p dx and C 2 = / E dx, 


which are also Casimir invariants of the Vlasov-Ampere equations. 
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Figure 3. Color map of S 4 (in A 4 units as defined in Sec. 12.311 given by Eq. (l33l) as a 
function of S 2 (in A 2 units) and S 3 (in A 3 units). 


3.2. Model with normal variables 


The second solution to Eqs. (fl9|) - (j20l) corresponds to the branch (/?, 0 , Rf R'f R'f ) = 
(1, 0, 2, 0) found in Sec. 12.31 and is given by 

R{Z) = i + e 2 - 

This leads to 

St = Si + 3 (33) 

*J2 

S5 = 2S 2 S 3 + 3 ( 34 ) 

These functions are plotted in Figs. [3]and|4j Unlike the model without normal variables, 
we can see on Fig. [3] that S 4 does not go to 0 as S 3 goes to 0. As a consequence, this 
model allows asymmetric distribution functions and hence seems to be more physically 
relevant. Furthermore, one can notice the difference in the amplitude of the closures 
between the two models (up to two orders of magnitude) by comparing Figs. |Tjand[ 2 ] 
and Figs. [3] and U The Hamiltonian and the Poisson bracket resulting from this closure 
are given respectively by Eqs. oa and (fill) with a and f3 given by Eq. (fT2l) and by 
replacing the closures S 4 and S 5 by Eqs. (1331) and (1341) . 

Concerning the global Casimir invariants, we show that this Poisson bracket 
possesses five Casimir invariants, as many Casimir invariants as field variables, given by 


Ci = I p dx, 

c 3 = / iX} 


C'2 — IE dx, 


C 4 = p 


S 3 


\/4S? + SI 


dx, 
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Figure 4. Color map of S 5 (in A 5 units as defined in Sec. 12.311 given by Eq. (l34l) as a 
function of S 2 (in A 2 units) and S 3 (in A 3 units). 


a = 


“ + fl' dx - 


From these expressions for the global Casimir invariants, we introduce the normal held 
variables p, M = u + pS 3 /(2S 2 ), Q 2 = p^/AS 3 + S 2 /S 2 , Q 3 = pS 3 /y/AS 3 + S 2 and E. 
Consequently, Bracket (fill) takes the particularly simple (normal) form 


{F, G} = 


Gm9xF p — Fm9 x G p + Att(GmFe 


FmG } 


2 G 3 d x F 2 - 2 GAF 


20 x r 3 


dx, 


Hamiltonian (flUll becomes 
1 


n = 


r2 usQa , p^2 


P M‘ -M^ + '-Qi + Ffj cb, 


and the Casimir invariants C 3 , C 4 and C 5 become 


C 3 = J Q 2 dx , C 4 = J Q 3 dx, C 5 = I M dx. 

As mentioned in Sec. ED one may be interested in using the pressure P and the heat 
flux q as dynamical variables, instead of S 2 and S 3 . Indeed, even though the reduced 
moments appear to be very convenient, their physical meaning may not be as clear as 
the usual pressure and heat flux quantities. The latter quantities can be expressed in 
terms of the reduced moments in the following way: 


p = p 3 S 2 = p 2 - 


P\ 


p ( 


0 


P C _ P 3 3 PiPx 2 | Pf 

q - 2 b3 ~ 2 2 P 0 + P 0 2 ’ 


in terms of which the closures take the form 
1 


Sa = 


P a 


P 2 Aq 2 


S* = 


Aq 
P 6 


P 2 q 2 

-b — 

p P 2 
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Expressed in terms of these variables, Bracket (fTIT) becomes 

3 P. 


{CG} = J 


G u d x F p — F u d x G p H- (G u d x F P — F u d x G P ) 

P 


+ 4 : 7 i(G u F E — F u G e ) H—( G U F P — F u G P )d x P -4 — -( G u d x F q — F u d x G q ) 

P P 

+ ~{G u F q — F u G q )d x q + p^OL2iF p G P + p°oi23F P G q + p 5 a32F q G P 
+ p 6 d33-F g Gg + p 2 P22G P d x (p 2 F P ) + p z ^2^G q d x {p 2 F P ) + p 2 P32G P d x (p 3 F q ) 


+ p 3 f3sG q d x (p 3 F q ) 


da;, 


where 


and 


a = d x 


4g/p 4 4g 2 /(p 5 P) - P 2 /(2p 6 ) 

6 q 2 /(p 5 P) ~ 3P 2 /( 2 p 6 ) 6 q 3 /(p 6 P 2 ) - 3Pq/p 7 


P = 


8 q/p 4 10 q 2 /{p 5 P) - 2 P 2 /p 6 

10g 2 /(p 5 P) - 2 P 2 /p 6 12 q 3 / (p 6 P 2 ) - 6 Pg/p 7 

Hamiltonian (TTUT) takes the simple form 


*4 


E 2 \ 

pu 2 + P + — j da;, 


and the equations of motion, obtained from c^P = {P, P}, are 


d t p = —pd x u - ud x p , 

dtu = —ud x u - d x P — P, 

P 

d t P = —ud x P — 3Pd x u — 2 d x q, 

d t q = —ud x q - 4qd x u - 2d x + j^d x (p 2 P 2 ) , 

d t E = inpu. 


We notice that these equations are identical (at least the ones concerning p, u, P and 
q) to the equations obtained with a bi-delta reduction [30, 31] [22, 32]. Therefore, as a 
by-product of our reduction procedure, we have proved here that the bi-delta reduction 
is Hamiltonian. This can also be shown by effecting a chain rule calculation relating the 
Vlasov-Poisson bracket pj9] to that of fluid streams [ 33 ] . 

A benefit of knowing the Hamiltonian structure of the reduced model is the ability 
to use the Poisson bracket to obtain the additional invariants, e.g., Casimir invariants, 
that can be tricky to derive directly from the equations of motion. For example, the 
global Casimir invariants C3, C 4 and C 5 for the present system are seen to be 

= / + dXi Ca = I da: ’ Cs = / (“ + f) da: ' 

We note that these invariants can be used to check the validity of numerical algorithms 
used for the integration of the equations of motion. 
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3.3. Comparison with Hamiltonian fluid models with 3+1 fields 


The same analysis done for 4+1 fields can be carried out for fluid models with 3+1 
fields, namely with the field variables (P 0 , Pi, P 2j E) or equivalently (p,u, S 2 , E). This 
was partly done in Ref. [29] (in the absence of electric field), where it was shown that 
Hamiltonian fluid models are given by closures S 3 that only depend 011 S 2 . This is also 
evident from Appendix C where the conditions given by Eq. (1B.7|) are automatically 
satisfied (since there is only one value for the indices). The fact that the closure for 3+1 
Hamiltonian fluid models only depends on S 2 is similar to the fact that the closures for 
4+1 fluid models are given by S4 and S 5 as functions of only S 2 and S3. 

A specific closure S 3 ( 62 ), which corresponds to the dimensional analysis performed 
in the present work, is given by 

c _ \ q3/2 

D3 — ad 2 , 

where A is a dimensionless constant. The Poisson bracket for this closure is 


{F,G } 3 = j \G u d x F f> - F u d x G, 


+ 4tt (G u F e -F u G e ) 


-(G u F 2 - F u G 2 )d x S 2 + 2AS 2 3/2 
P 


G, 


d x 


dx. 


A 

p \p J p V p 

It should be noted that the dimensional analysis provides a family of models (labeled 
by A). However there are only three fundamentally different models: one for A = 0 and 
the others for A = ±1, since all of the other models can be rescaled to A = ±1 by an 
appropriate rescaling of S 2 , e.g., S 2 = S 2 / A 2 . Moreover, the two models A = 1 and 
A = — 1 are linked by symmetry [29j. The model for A = 0 has the two following global 
Casimir invariants: 


Ci = p dx 


and 


C 2 = E dx, 


in addition to the family of Casimir invariants 


C= pn(S 2 ) dx, 


for any scalar function k. The two Casimir invariants Ci and C 2 are identical to the ones 
for 4+1 fields. Concerning the model with A = 1, the Poisson bracket with 3+1 fields 
has Ci and C 2 as Casimir invariants, and also has two additional Casimir invariants 

C 3 = j pS\^dx and C 4 = J (u — 2 pS^ 2 ^ dx. 

Therefore, in total it has four Casimir invariants, i.e., as many as the number of field 
variables. The common point between this 3+1 model with A = 1 and the 4+1 field 
model with normal field variables is that both have a generalized velocity as Casimir 
invariant. It should also be noticed that the 3+1 fluid model has one Casimir invariant 
of the entropy type, i.e., of the form J p<f>(S 2 ,..., S'jv) dx, whereas the 4+1 fluid model 
has two of this type. 
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4. Summary 

In summary, starting from the one-dimensional Vlasov-Ampere equations, we built two 
Hamiltonian models with the first four moments of the Vlasov distribution function and 
the electric field as dynamical variables. Our reduction method relied on the preservation 
of the Hamiltonian structure of the Vlasov-Ampere model. The closures we obtain 
were derived from a dimensional analysis argument. We showed that there are only 
two Hamiltonian closures obtained by this method. A fundamental difference between 
these two models was characterized by their Casimir invariants: one model has only 
two global Casimir invariants (preserved from the Vlasov-Ampere system), whereas the 
second model has three additional ones, two of the entropy-type and one generalized 
velocity. 
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Appendix A. Independence of a and f of Bracket (llljl on p, u, and E 


In this appendix, we consider the following bracket defined on functionals of the form 
F\p,u,S 2 ,...,S N ,E] for N >2: 


{F, G} = 


G u d x F p - F u d x G p + 4 n(G u F E - F U G E ) - ~(G u F t - F u Gi)d x Si 


F l G j n FA n Gj , 

+ otij -- + d x — fij— da, 

P P \pj P 


(A.l) 


where a and f are matrices satisfying f t = f and d x J3 = a + cd, assuring antisymmetry 
of the bracket. Here we assume a priori that a and f depend on both the dynamical 
variables p, u, Sk (for k > 2) and E, and their derivatives dfp, d™u, dfSk and dfE for 
n > 1. Repeated indices are implicitly summed from 2 to N , unless specified. We seek 
necessary conditions on a and ft for Bracket (lA.lj) to satisfy the Jacobi identity, 


{F, {G, H}} + {H, {F, G}} + {G, {H, F}} = 0. 


In this appendix, we prove that a and ft do not depend on the variables p, u and E and 
their derivatives dfp, dfu and df E for n> 1. 

First we split Bracket (1A.1D into two parts 


{F,G} = {F,Gy + {F,G}', 
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where the first part, 


{F,G} J = 


G u d x F p - F u d x G p + 4t t(G u F e - F U G E ) - ~(G u Fi - F u Gi)d x Si 

P 


dx, 


satisfies the Jacobi identity puns]. The Jacobi identity is then equivalent to 

{F, {G, H} J Y + {F, {G, H}*} J + {F, {G, HYY+ 0 { f,g,h)= 0, (A.2) 

where C)(f,g,h ) denotes the summation over circular permutation of any three functionals 
F, G and H. Using the lemma stating that only the functional derivatives with respect 
to the explicit dependence on the variables need be taken into account for the Jacobi 
identity fH], the first term becomes 


{F,{G,H} J Y = 


&ij Fi 

L P P 

+^d, 

p 


Gj 

P 


H u d x ( — - G u d x — + —d x H u - L d x G v 


Gj 


3 d x H u - 


P 


H± 

P 


Hj 

P 


d x G u 


Gj 


Hj 


P 

dx, 


P 


(A.3) 


where we have used the fact that f is symmetric. The second term in Eq. (1A.2I) is 
{F, {G, «}’}■’= f [{G, HY&F r + {G, H}*d x F u + 4*({G, HYJY - {G, H}* E Fj 


- -YC.HYJj- F,{C,HY\d,S, 


where 


{G,HV = 


a 




Ot; 


P P 

+ (-ir<% 


Gj Hj fijj 

P P PI 


0a j,j Gj 


+ 


d x 

dfiij 


dx, 


57) — - a 

P J P 
G 


(A.4) 


H± 

P 


d. 


dd'fp p dd^p V P 


P 


{G,HY u =(-l) n d. 
{G,HY k =(-l) n d. 
{G,HY E = (-l) n d : 


datj Gi d/3. 


dd™u p dd'fu 
dfiij 


"41- 

P 


Hj 

P 


da jj Gj 


+ 


ddfS k p dd'fSk 
dfijj 


a,|U 

P 


Hj 

P 


daij Gi 


+ 




Hj 


dd x E p dd'fE \ p j j p 
We consider the terms of the type (F U1 Gj , Hj) in the Jacobi identity f1A.2[) . These terms 
only come from Eqs. (IA.3[) and (1A.4I) . By using successive integrations by parts and 
assuming that the boundary conditions are such that the associated boundary integrals 
vanish, the Jacobi identity for these terms becomes 


da.jj G j 


+ 


df 


v 


-OLjjdx 


dd'fp p dd'fp 
Gj\ Hj F u 

->T~ 

dfiij 


d x 


Gi 

P 


^F X - ^±d x 


P 


P 


Gj F„ 


OLi 


P P 


47r 


daij G 


i] dfjj 


+ 


dd'fE p ddfE 


dr 


Gj_ 

P 


H 


1 d n F 

U x r u 


p 


(A.5) 


da.jj G 


+ 


8 , |* 


88 jS t p 8 &jS k V P 


H 


P 


P 


—d x — d x S k > dx+ 0( F ,g,h)= 0. 
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Choosing F = f u dr, G = f pSi dr and H = pS m , Eq. (IA.5D leads to the necessary 
condition 


1 rj &lm rj . dotlm p,n f 1 r, r, 

-o x aim = —d x p + d x -d x S k 
P P- Od^S k \p 


(A.6) 


However, we have by definition 

8*<*, m = ^ep'p + ^dr'u + 


da 


lm 


ddfp 


ddjtu 


o: +1 s k + 


da 


lm 


dd?s k x dd?E 


df +1 E + 


da 


lm 


dx 


where the summation over n is implicit and dai m /dx denotes the derivative of cq m with 
respect to its explicit dependence on x. Eventually, Eq. (IA.6[) writes 


daim orj i | dai m 


P \dd%p 

&lm 




da, 


lm 


dd!ts k 


d n x +1 S k + 


da 


lm 


dd™E 


df +1 E + 


da 


lm 


dx 


2 d *P + 
p 2 dd™S k \p 


(A.7) 


By canceling the only term that depend on df +1 p in Eq. (1A.7I) . we can show that 
dai m 


dd x p 


= 0. 


By performing an induction on v down to v = 0, we can show that a does not depend on 
p and its derivatives. Because the dynamical variables are independent, using the same 
reasoning we prove that a cannot depend on v, E and their derivatives, nor can it depend 
explicitly on x. The same result can be obtained for (3 by choosing G = f pSix da;. 
Therefore a necessary (however not sufficient) condition for Bracket (1A.1I) to satisfy the 
Jacobi identity is that a and fl do not depend explicitly on x , p, u and E, as well as the 
derivatives &£p, d x u and d x E for all n G N. 


Appendix B. Dependence of a and fl of Bracket (lllll on S k 

In this appendix, we derive some necessary conditions on the dependence of a and fl 


(and their derivatives) of Bracket (fITjl on S k . Following Appendix A, we consider two 
sets of functionals 


(F, G,H) = [ ux dx, / pSi da;, pS n 


and 


(F, G, H) = ux dx, / pSix dx, pSn 


which we insert into Eq. (1A.5I) . Thus we find the necessary conditions 

dfllm rvn 


daim an 

a im = n 0 d x S k , 


n- 


-d x S k = 0, 


(B.l) 


ddSS k x ddfS k 

for all l and m, where we recall the implicit summation over repeated indices. We 
assume that a and fl depend on the derivatives of S k up to order u, where 


v = max{n G N s.t. da/dd^S 0 or dfljdd n x S 0}. 
































Hamiltonian closures for fluid models with four moments by dimensional analysis 19 


From the first of Eqs. (IB . 1 1) we have 


dxPim — ^ ^ 


dfllm 
n ddfSk 

n =0 


d' +1 S k = ai m + 


E" 


&ml s ^ 
72—0 


9c^lm ^ dot-mi 


dd™S k dd!tS k 


d:s k . 


(B.2) 


Differentiating Eq. (IB. 21) with respect to d v+1 S 3 leads to 
dflir 


dd^Sj 


= 0 . 


As a consequence, the highest derivatives of S appear in a; thus v becomes 

v = max{n e N s.t. da/dd x S 0}. 

The Jacobi identity (IA.2|) reduces to: 


{F,{G,H}}+0 { f,g,h) = 


„ F t {G,H}* f o {G,H}* 

Qtij + Pij 

. 9 


Or. I - 
P 


dx 


(B.3) 


P P 

+ C)(F,G,H)= 0. 

This identity corresponds to the Jacobi identity for the subalgebra of observables 
F\p, S 2 , ■ ■ ., Sn]. Expanding Eq. (IB.3|) gives 

an Fi 


{F,{G,H}}+ O (f , g ,h)= 


P P 


da k i G_k dflki q (Gk 

dd™ Sj p ddf. Sj x V P 


p 


+ (-i) n+1 -d x 

P L P 


Pij Qn 


dotki G k , dflki 


+ 


+ d " 


L dd^Sj p dd™Sj 

a ij Hi | Pij q f Hi M d(Xi k Fi G k 


Gi 


dr — 


p J J p 


H, 


P P 


P 


P 


ddfSj p p 


-dr ( d 
P 


a lj H l + PfiQ (Pi 


+ d: 


L P P 
oikj G k /3fc 


P 


P 


+ i A 


L P P 

Choosing consecutively 


P 


P 


dfiik Gk 
dd x Sj p 
da u Hi dpi 


+ 


ddfS 3 p dd^Sj 

(■ F,G,H) = I pSi , / pSk dx, / pSi dx^) , 


dr I — 

P 


Fi 

P 


— \ dx. 


(F, G,H)= pSi, / pS k x dx, / pSi dx , 


and 


(F, G, H) — ( pSi, I pS k dx, / pSvx dx ) , 

dan 


we get the following three conditions: 

doiik dfiik 


dd x Sj dd^Sj 
„ dan 

^ kj ~ddpS 3 = °’ 
dflik 

Plj dd" '.S' 


OLlj + ^kj 


dd"S 3 


= 0 , 


(B.4) 

(B.5) 

(B.6) 
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Due to the fact that d x (3 = a + a *, Eqs. (IB.51) and (IB.61) are redundant. We assume that 
a depends only linearly on df.Sk- We show in Appendix C that this is the case for fluid 
brackets. As a consequence, we write 


aij(S, d x S, ...,9fS) = A tj (S , d x S ,..., df-'S) A 7 d x S ,..., df~ 1 S)dfS, 


By inserting this expression into Eq. (IB.41) . for the Jacobi identity we need to impose 


r iljm 


dfiik 

likj ddfr'Sj 


A 7 kjm'Ylij 


0 , 


for all (i, k, l, m) to make the term proportional to dfS m vanish. However, thanks to 
Eq. (fR~2l) we have 

dfik 

ddf- l S, likj ~ ~ :kir 

This eventually leads to the following conditions: 

7 ijmT/kij Ikjmllij- (B. 7) 

These commutation relations remind us of the conditions for Lie-Poisson brackets based 
on Lie algebra extensions to satisfy the Jacobi identity of Ref. [33] • These conditions on 
the tensor 7 are necessary but not sufficient. 


Appendix C. Jacobi identity for fluid models 


In this appendix we fold necessary and sufficient conditions for the Jacobi identity 
for Bracket OB. We start from the one-dimensional Vlasov-Ampere bracket given by 
Eq. (|S]) and perform a change of variables, from / to (p, u, S n > 2 ) defined by 

p= I fdv, pu= j vf dv, p n+1 S n = j (■ v-u) n fdv. 

Using the following chain rule expression for the functional derivative with respect to /, 


v — u - 

Ff — F p 4 - F u + F n 

p 



n + 1 S n _ 1 (v - u) 

-— r> - 


and after some algebra, we show that the Poisson bracket 
a and f3 given by 


reduces to Eq. (IA.1|) with 


c^nm n() x S r) . m 1 nS n -id x S m ? 7 .(m V n7nS m —id x S n , (C.l) 

finm {jYl ~\~ 7lfS n - )-m —1 TYl{jl A 1) S n S m — | n(jn A l)S m S n -l, (C.2) 

where n,m> 2 and S\ = 0. The resulting bracket is of Poisson type. Next, we truncate 
the matrices a and j3 such that a mn = 0 and /3 mn = 0 for m > N or n > N. The matrices 
a and f3 depend on S n for n — 2 ,..., 2 N — 1 . We restrict ourselves to the case where a 
and f are functions of (S 2 , • • •, Sn), i-e., we introduce N — 1 closures Sk = Sk(S 2 , ..., Sn) 
for k — N A 1, ■ ■ ■ ,2 N — 1. In this truncation/reduction, the bracket is no longer of 
Poisson type in general. In this appendix, we establish the necessary and sufficient 
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conditions for the Jacobi identity to be satisfied. From Appendix A and Appendix B 
this Jacobi identity is seen to be 


{F\ {G, H}}+ O (F,G,H)~ 

dan 


— dk 


3 dd x Sj 


+ 


0)j d x 


da 


kl 


dd x Sj 


d x . 

P J P 

m Gkm 

p p p p L 

otijd x 


-d x 


a 


v 


Gk 

P 

da H 

dS, 


p 


o dfki q da k i 

Pij a c A? 


dSj 


dd x Sj 


da ik da h 

+ a,i ~as~ + ati as~ 


da 


ik 


dd x Sj 


® k jd x 


da,. 


dd x Sj 


Fi 1 


d x 


P P 
fkjd x 


ga h 

p ) p 
da u 


dfki 


dai 


aij ^ ^ dk i ^ ~ 


da 


kl 


dS , 


J dSj 


v dd x Sj 


a,j 


da, 


ik 


dd x Sj 


dd x Sj 


+ d x 


F, \ 1 G, II, 


dij d x 


da 


kl 


dd x Sj 


a,j 


P J P P P 
da ik dan 


d 


da 


kl 




dSi 


+ Ol,j 


ddik 
dSi 


dd x Sj 


a k j 


dd x Sj 


H— d x ( 
P 

- dijd x 


(HA F\Gk_ 
\P ) P P 
daik 
dd x Sj 


A. 


da 


ik 


. J dSj 
dan 
ak3 dd x Sj 


da k , ddii 

aij dd x Sj +akj ~dS~ 


+ d x 


FA 1 


-d x 


P J P 


El 

p 


q ddik q da k , dai k 

p lj ~QS~ ~ P ij ~Fa~a A? 


dd x Sj 


dd x Sj j 


P 


G, 


P 


+ -d x [ — \d x \ — \- 


II, 


P J P 


dk_ 


ddii 

dSi 


A. 


da 


ik 


3 dd x Sj 


dk 


Gk 

P 

dan 


3 dd x Sj i 


(C.3) 

dx. 


Choosing F = pSi, G = f pS k dx and H = f pS, dx, Eq. (1C.3|) reduces to 


a^j / 


da 


Ik 


p V dSj 


- d x 


da 


Ik 


dd x Sj 


+ d x 


Pij 

P 


d x 


da 


Ik 


dd x Sj 


da 


Ik 


dS, 


O kj ddii 

—~dS~ 


akj dan . f dan ^ . ®-ij da k i da k i 

H- -^77- + -7T7 r-Z~d x a k j + -7777- + 


3*( y 1=0. 


p dSj 1 pdd x Sj~ J r> ' p dS, 1 Hr),3, J \ p J ~ (C ) 

Equation (1C.41) can be split into two terms with only one depending on p. To make the 
term that depends on p vanish, we have to impose 


A 




d x 


da 


Ik 


dd x Sj 


da 


Ik 


dS. 


3 


ddu , da a . da ki 

+ Ad — c _Qi b - °> ( C - 5 ) 


dSj 


dd x Sj 


dd x Sj 


for all (i,l,k). In addition, canceling the term in Eq. (1C.41) that does not depend on p 
leads to 


a. 


13 


da, k „ da, k 
- d x 


dSj 


^ + a kj f 


da a r, da a 

- o x - 

dSj 


dd x Sj 


Fan 


dd x Sj j 

for all (i,l,k)- With these constraints, Eq. (1C.31) becomes 


da ki da k i 

- a .. 

dSj ‘ 


dd x Sj 


= 0, (C. 6 ) 


{F,{G,i/}}+ 0 (F ,G,H) = 


FA 1 


d x - \-d x — — 


P J P 


GA Hi 


P J P 

-A 


A 


V 


ddki 

dS~ 


da k i da h 

— Pk 


13 dd x Sj 


3 dd x Sj 
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+ d x 


FA 1 


-d x 


P J P 
G k 


H,\ Gi 


P 


P 


-dA — \dA — \- 


p 

Hi 


P L 


A. 


dfik 

dS, 


- fi- 


da 


kl 


- Aj 


dai 


ik 


p / p 


13 dd x Sj yj dd x Sj 

dfu da ik dan 

dG dS ■ A? A A C dkj 


dd x Sj 


dd x Sj j 


(C.7) 

dx. 


Choosing F = pS il G — f pS k dx and H = f pSix dx leads to 


-d x 


1 (_ dB ik n da k i _ da ik 

Pij^TT ~ ~ A. 


LP 








= 0 , 


which has to be satished for any p, and therefore 


df ik da ki 

Plj r } 0 A 


0 S,- 






(C. 8 ) 


for all (z, /, P). With this additional constraint, Eq. flC.7j) is always satisfied, which proves 
that Eqs. (1C. 511 . (1C.61) . and (1C. 8 ft are necessary and sufficient conditions for Bracket (HID 
to satisfy the Jacobi identity. 

By introducing the expressions of a and ft given by Eqs. (1C. 1 jl and (1C.2D into 
Eqs. (1C.5D . (1C.61) . and (1C. 81) . we end up with the following constraints : 

(C.9) 
(C.10) 


P _ P 

1 iklm L ilkmi 

^ikl ^Ikii 


where 


_ rA 

iklm. hW 


(1 — i — l)Si + i_ 2 + jSj-i 


dS l+ i\ 


dSi 


rfc—1 


(i + l)S i+l . 1 -(j + l)S j 


dSi + i -1 


dS i+ i_ 1 SS'fc+j-i 


and 


AC 

^ki = [(/ + j) A+i -1 - j(l + l)SiSj - 1 - l(j + l)^-5,_i] 

+ l(i + P)S , /_i»5’j_|_fc_i + (/ + l)(z + A; — l)S , ;S , j_|_fc_ 2 - 
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